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ABSTRACT 

Amplitude  and  phase  relationships  between  sinusoidal  pressure  dis- 
turbances in  the  exit  plane  of  a  conical  nozzle  and  mass  flux  response 
were  experimentally  determined  for  the  subsonic  flow  of  air  in  the 
nozzle.  The  experimental  results  were  compared  with  values  predicted 
by  using  the  modified  asymptotic  WKB  method. 

Experimental  measurements  were  made  with  a  constant  temperature 
hot  wire  anemometer  and  a  condenser  microphone.  Disturbance  frequen- 
cies from  20  to  250  Hz  were  investigated  for  Mach  numbers  up  to  .34. 

The  experimental  values  of  amplitude  factor  agree  quite  well  with 
analytical  predictions  for  all  Mach  numbers  tested.  Experimental 
values  for  phase  angle  differed  from  those  predicted  but  exhibited  a 
si  mi  iar  v»renu . 
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NOMENCLATURE 


a 

speed  ot  sound 

ao 

speed  of  sound  at  stagnation  condit 

A 

flow  area 

A* 

flow  area  where  M  =  1 

cr 

C2 

arbitrary  constants 

D 
Dt 

substantial  derivative 

F 

defined  in  Eq.   (24) 

9 

defined  in  Eq.   (29) 

h 

enthal py 

H 

stagnation  enthalpy 

i 

yfT 

M 

Mach  number 

H,. 

N2,   N3,   N4 

defined  in  Eq.   (62) 

N 

amplitude  factor 

P 

Pressure 

q 

velocity  (magnitude) 

q 

velocity  vector 

r 

radius  measured  from  apex  of  cone 

R 

radius  where  M  =  1 

s 

entropy 

S*' 

modified  Strouhal   number  oor/q 

T 

temperature 

To 

stagnation  temperature 

t 

time 

UT 

Vl 

real  and  imaginary  components  of  Z, 

u«»  v2  real  and  imaginary  components  of  Z2 

V  hot  wire  voltage 

V  hot  wire  zero  voltage 
o 

w 

X 

y 
zrz2 

a 
6 
Y 
C 

n 
u 

IT 
P 

t  dimension! ess  time  a  t/R 

<J>  dimension! ess  enthalpy  perturbation  H  /a 


$ 


frequency  number  wR/aQ 

=  M2 

defined  in  Eq. 

(65) 

defined  in  Eq. 

(34)   * 

q/q0 

a/aQ 

ratio  of  speci 

fie  heats  (= 

1.4  for 

air) 

r/R 

defined  in  Eq. 

(25) 

(pq)'/pq 

• 

transformed  5 

p'/p 

density 

dx 
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^  phase  angle 


oj  angular  frequency 
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■   I.  INTRODUCTION 

The  engineering  application  of  nozzles  in  the  control  of  fluid  flow 
is  extremely  widespread.  It  ranges  in  scope  from  small  fluidic  devices 
to  large  rocket  thrusters.  In  the  steady  state  there  exists,  for  a 
given  absolute  pressure  ratio  across  the  nozzle,  a  corresponding  mass 
flow  rate.  The  introduction  of  pressure  fluctuations  at  the  nozzle 
boundaries  affects  the  flow  and  certain  dynamic  characteristics  are  ex- 
hibited for  a  given  set  of  conditions,  i.e.,  nozzle  geometry,  type  of 
fluid  medium,  and  the  nature  of  the  pressure  fluctuations.  Because 
Such  fluctuations  are  quite  common,  especially  in  systems  possessing 
cyclic  operating  components,  a  knowledge  of  the  dynamic  characteristics 
becomes  important.  In  the  past,  a  quasi-steady  relationship  between 
mass  flow  rate  and  pressure  fluctuations  was  assumed.  Thus,  an  unsteady 
situation  was  reduced  to  a  sequence  of  instantaneous  steady  state  cor- 
respondences. The  validity  of  this  assumption  is  questionable. 

Oppenheim  and  Chilton  [3]  conducted  a  literature  survey  on  pulsating 
flow  measurement.  Most  of  the  investigations  noted  had  to  do  with  the 
effect  of  pressure  fluctuations  on  the  time  average  of  flow.  Elrod  [1], 
working  in  connection  with  a  stability  analysis  of  externally  pressurized 
thrust  bearings,  performed  a  frequency  response  analysis  to  predict  the 
amplitude  and  phase  relations  between  the  instantaneous  flow  rate  and  a 
sinusoidally  fluctuating  exit  pressure  in  an  infinite  conical  nozzle. 
Although  the  chosen  geometry  was  extremely  simple,  Elrod  [1]  stated 
"It  is  believed  that  the  essential  characteristics  of  a  compressible 
fluid  expanding  and  accelerating  from  a  stagnation  reservoir  are  well 
enough  realized  for  the  results  to  have  semi -quantitative  practical 
significance.."  He  further  assumed  that  the  flow  was  adiabatic, 
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friction! ess,  and  irrotational.     The  analysis  using  the  WKB  method  is 
asymptotically  valid  for  high  frequency  number  W=coR/a  and  low  Mach 
number.     Unfortunately,  the  numerical  results  were  wrong  because  of 
a  wrong  choice  of  sign  for  the  divergence  term  of  the  developing  equa- 
tions. 

Chiang,  et  al.,  [2],  realizing  the  generality  of  Elrod's  analysis 
(it  can  be  applied  to  either  converging  or  diverging  flow),  corrected 
the  choice  of  sign  and  extended  the  region  of  validity  to  lower  fre- 
quency and  higher  Mach  number  by  modifying  the  WKB  method.     They  also 
investigated  the  problem  of  pulsating  flow  in  a  finite  conical   nozzle. 

The  objective  of  this  thesis  was  to  experimentally  investigate  the 
mass  flux  -  pressure  relationship  for  pulsating  subsonic  flow  in  a 
conical   nozzle  and  to  compare  experimental  results  with  values  pre- 
dicted using  the  modified  WKB  solution. 

II.     THEORETICAL  ANALYSIS 

Elrod  [1],  in  his  original  analysis,  assumed  an  adiabatic,  friction 
free  (i.e.,  isentropic)  flow  originating  in  a  stagnation  reservoir. 
Under  these  conditions,  the  flow  remains  irrotational   in  accordance 
with  Kelvin's  theorem  and  the  following  governing  equations  pertain: 

Momentum 

b&  +±.  vf  ^O  0) 

D±       ? 

Continuity 

J2p  +pV-q.  =  0  (2) 


which  can  also  be  written  as 


Dt1  °  l  ' 


Gibb's  Equation  for  a  simple  compressible  substance 


Tc/s  ^  c/h--i-  d-p  =  O  (isentropic)        (4) 


For  irrotational  flow: 


(<VV)%  afW$'%)  =  V-f-  (5) 


2 

(see  Robertson  [4],  p.  46.) 

Using  Eq.   (5),  Eq.   (1)  can  be  rewritten  as 


|i  =  -|Vf-V-^  (6) 


Stagnation  enthalpy  is  defined  as 


U*Jx+-3L  (7) 

z 

Equations  (4),   (6),  and   (7)  combine  to  yield 


-ff-'-VH  (8) 

The  scalar  product  of  q"with  Eq.    (1),  together  with  Eq.   (4),  gives 
2JL   --JLqL-VP^-qsVh  =  -M-  +  -2A.  ,  (9) 
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which  can  also  be  written  as 


_DH  ,  J_  d£. 

by  considering  a  time  derivative  form  of  Eq.  (4) 
The  acoustic  velocity  is  defined  by 


a  = 


mi 


Therefore,  by  the  chain  rule,  Eq.   (10)  becomes 
±  JUL  -   -L  .££.  -    tJUo 


00) 


OD 


02) 


Taking  the  total  derivative  of  Eq.   (12)  gives 


03) 


but 


D      ZJUp  .    B     d>Up  +^(&-V^p\-d%z-*V^  P 


=^te^)-H'^c 


04) 


Substituting  Eqs.   (3),   (8),  and  (14)  "into  Eq.   (13)  above  produces 
El  rod's  [1]   "propagation  equation"  for  stagnation  enthalpy: 


1  DH_ 


=  \raH+VH-7^p 


05) 


For  spherically  symmetric  flow  in  a  conical  passage  with  slight 
flow  deviations  which  are  also  spherically  symmetric,  the  flow  is  des- 
cribed by  only  one  spatial  coordinate  -  namely  the  radius  r  originating 


11 


at  the  apex  of  the  cone. 

Chiang,  et  al.,  [2]  used  linear  perturbation  theory  to  express  a 
quantity  as  the  sum  of  a  steady-state  value  and  a  perturbation. 

Stagnation  enthalpy  was  written  as 

H  =  H  +  H'  (16) 

where  the  prime  denotes  a  perturbation  term. 

Since  the  stagnation  enthalpy  is  uniform  throughout  the  mean  flow, 
it  is  apparent  that 

VW-0  07) 

and 


DH 


*   O 


Hence,  Eq.   (15)  may  be  written,  to  the  first  order  in  small   quan- 
tities, as 


=  V^'-rVH'-V^f 


As  mentioned  previously,  spherically  symmetrical  flow  is  one- 
dimensional  .  Chiang,  et  al.,  [2]  pointed  out  that  for  this  case 

dr 


08) 


09) 


(20) 


The  upper  sign  pertains  to  diverging  flow  and  the  lower  sign  is 
related  to  converging  flow.  This  choice  of  signs  is  based  upon  the 
fact  that  for  diverging  flow  the  velocity  and  increasing  radius  are 
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identical  in  direction.  For  converging  flow,  the  directions  are  opposed, 
Therefore,  the  lower  (negative)  sign  applies  to  nozzle  flows. 

For  convenience,  El  rod  [1]  introduced  the  following  dimensionless 
variables: 


<t>- 


_  H 


) 


X  -  r 
5"  R 


(21) 


where  a  is  the  acoustic  velocity  at  stagnation  conditions  and  R  is 
the  fictitious  radius  at  which  M  =  1  for  conical  flow.  In  terms  of 
these  new  variables,  the  "propagation  equation"  Eq.  (19)  becomes 


dt[ 

\a  J  t>r 

a0 

where 

D 

DT 

-3      +  . 

ar 

4-  i_ 

CL0    ^ 

la 


(22) 


(23) 


Elrod  [1]  considered  sinusoidal  disturbances  and  found  it  bene- 
ficial to  perform  a  frequency  analysis  by  representing  $  in  the  form: 


^  =  F0Oe 


iwr 


(24) 


W  is  a  dimensionless  frequency  number. 
It  was  also  assumed  that 


Fft)«d-)*^{txw/M<i5}V!) 

The  resultant  differential  equation  for  n(S)  is 


(25) 


(26) 
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where 


=  -&. 


a. 


P  = 


-    a0 


(27) 


?-/i&d* 


(28) 


and 


m  -  -  £  *.&■)*  ♦&  /-ft)1 


d? 


-    M' 


OWL 


2  Z  Z 


M' 


(29) 


Using  the  steady  state,  one-dimensional,  isentropic  relationships 
for  area  and  temperature  as  functions  of  Mach  number  (see  Shapiro  [5]), 
the  following  relationships  hold 


p 


•*  W.  • 


I  +£1^ 


C30) 


A 


2 


=  h-^  o*^y 


r+i 


4(y-0 


(3D 


-  Chiang  [2]  arrived  at  the  general  solution  for  Eq.  (18)  by  using 
the  WKB  method.  (See  Morse  and  Feshback  [6]). 
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(32) 


^(5)»C1e*f[-iw!+-^jyi)df] 


C,  and  C«  are  arbitrary  constants.  Eq.  (24)  now  is  written  as 

f=(ff  e^[t x. w/MdfJ  j"  c,  c^ [iw(r-f )] 

Chiang  [2]  states  that  the  C,   term  with  the  exp  ||  iw(x-£)J     indi- 
cates an  outgoing  wave  (in  the  positive  r-direction)  while  the  Cp  term 
with  exp    Liw(T+£)~]    indicates  an  incoming  wave;  the  term  exp      iwf  Md^  J 
represents  a  mean  convection;  and  the  factors  exp      +  -*-  J    g(£)d£  I 
represent  attenuation  of  waves.     For  the  case  of  an  infinite  conical 
nozzle,  with  pressure  fluctuations  impressed  at  the  throat,  C2  =  0. 

The  WKB  solution  is  asymptotically  valid  in  the  limit  of  small  Mach 

3  2 

number  M  and  high  frequency  number  W  with  an  error  of  the  order  (M  /w  ). 

Chiang  [2]  extended  the  solution  to  larger  Mach  number  and  smaller 
frequency  number  W.  A  summary  of  this  extension  follows: 


(33) 


15 


Assume 


(35) 


^=C|eiwS2,(5)+caeiw?zaC?)  (34, 

After  substituting  into  the  previous  differential  equations  for  n 

and  making  comparisons,  asymptotic  approximations  of  Z-,  (£  )  and 

Z0(E  )  are  determined. 
do 

where  £     is  sufficiently  large  (i.e.,  far  enough  away  from  the  throat) 
so  that  the  WKB  method  holds.     Differentiating  Eqs.   (35)  and  (36) 
yields 


2,'feJ-aB  &-) 


(37) 


Kd)-ix  fit.) 


(38) 


It  is  further  shown  that  Z-,  and  Z?  are  complex  conjugates  and  that 
Eqs.  (35)  -  (38)  can  be  used  for  initial  conditions  to  solve  the 
differential  equations  for  Z-,  and  Z«. 
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Let 


Z,=   U,+X\^ 


Z,=   U^+Avj 


then 


(39) 


U,=   U^ 


"Uf«\5 


(40) 


After  changing  variables  by  using 


x=  m 


and 


$  = 


-    d* 


*n 


Chiang  showed  that 


$=-hi 


and  arrived  at  the  following  differential   equations 


*?(*)+■£ 


-.gM  u,W  =0 


V'(*)Hrjj|f  VM-2wu;(x)j 


(41) 


(42) 


(43) 


(44) 


(45) 
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with  the  initial  conditions 


u;(^  =  o 


(46) 


where  x  is  the  value  of  x  corresponding  to  E,   =  £  .  3,  C,  *,  and  g 
have  all  been  previously  defined  as  functions  of  M  and  are  therefore 
functions  of  x. 

It  is  always  possible  to  select  a  Mach  number  M  small  enough  to 


satisfy  the  restriction  on  the  WKB  solution, 

(47) 


W2* 


The  amplitude  and  phase  relations  between  pressure  and  mass 
flux  perturbations  are  determined  in  the  following  manner.     Eqs.   (8) 
and  (10)  are  used  to  obtain 

1%L=    T    JH.'  (48) 


JL  ^Y  -    Zj£  ±    9v-^i^  (49) 

f     t*     ""      die  ^r 
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The  mass  flux  perturbation  is  written 

M'  =  e*/  +  %e'  (50) 

After  introducing  the  dimension! ess  pressure  and  mass  flux  perturba- 
tions: 


ir=^-  k  (si) 


Chiang  [2]  showed  that 

BIT  DJU&  +  .M  ,     $J*-£. 

Furthermore,  since  both  n  and  y  are  periodic  in  time  and  contain  the 
factor  exp  £iwt] 


2L 


9TT 

(54) 


For  the  case  of  an  infinite  conical  nozzle  with  sinusoidal 
pressure  fluctuations  impressed  at  the  throat,  C~  in  Eq.  (34)  be- 
comes zero  and 
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^?)=C,e^[-IW£]  (u1+xur) 


therefore 


*•$$)***+[-*-»  fn<tf] 


a^ 


3? 


i.-/$/LL-JfciJ\  +  JE (u  ^3 


where 


=  phase  angle  =    arctan  /  ^L  ^3      fy  % 


(55) 


L  J  (56) 


x  e<op[^w(r-  j)l  (U|  +xv;) 

Taking  the  logarithm  of  $  and  differentiating  gives 


^i=iW  (57) 


(58) 


Substitution  of  Eq.   (57)  and   (58)  Into  Eq.   (53),  together  with. 
Eq.   (54)  yields 

3L  =  is  JTMljve^  C59) 


/V  =  amplitude  factor  =/~r 7T  1  (50) 


(61) 


20 


and 


(62) 


Ws=-A/, 


1V+  =  -  Wa  -  W  x"  * 


A  dimension! ess  parameter  designated  as  the  modified  Strouhal 
number  S*  was  defined  by  Chiang  [2] 

<-*  _  60  r 

%  C63) 

and  shown  to  be  related  to  the  frequency  number 


W  a 


-   M  S*  (64) 


I3  ? 


For  a  given  Mach  number  M,  it  is  more  convenient  to  express  the 
mass  flux-pressure  relationships  in  terms  of  the  modified  Strouhal 
number  S*  rather  than  the  frequency  number  W. 
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(65) 


III.     NUMERICAL  SOLUTION 

Equations  (44)  and  (45)  are  transformed  to  a  system  of  four  first 
order  differential  equations  by  defining  the  variables: 

y^x)  =  u^x) 

y2(x)  =  u{(x) 

y3(x)  =  v^x) 

y$M  -  v{(x) 

The  resultant  system  of  differential  equations  is 


(66) 


*—*&  +  '*"*]+ sp* 


y4 


A  numerical   solution  for  this  system  of  equations  was  obtained 
using  Hamming's  Modified  Predictor-Corrector  Method  of  numerical   inte- 
gration.    (See  Ralston  and  Wilfe  [7]).     It  is  a  stable,  fourth  order 
integration  technique  which  uses  four  preceding  points  for  the  evalu- 
ation of  each  successive  point.     PROGRAM  TWO  makes  use  of  an  IBM  [8] 
subroutine,  DHPCG,  to  perform  the  numerical   integration. 

Evaluation  of  the  initial  conditions  v-,  (x  )  and  v,  (x  )  is  accom-  ■ 
pi i shed  using  PROGRAM  ONE.     Subroutine  GOPHI  evaluates  the  function 
g(x)/$.     Numerical   integration  of  v. (x  )  is  performed  using  a  twelve 
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point  Gaussian  Quadrature  formula  which  integrates  polynomials  up  to 
degree  23  exactly.  Subroutine  DQG12  is  used  for  this  purpose.  (See 
IBM  [8],  p.  299.) 

With  u,,  u, ,  v,,  and  v,  known,  subroutine  OUTP  of  PROGRAM  TWO 
computes  values  of  amplitude  Factor  N  and  phase  angle  ip  for  a  given 
Mach  number  M  and  frequency  number  W.  The  corresponding  modified 
Strouhal  number  S*  is  also  computed. 

The  programs  have  been  written  for  the  case  of  air  with  y  =  1.4. 

IV.  EXPERIMENTATION 

A.  EXPERIMENTAL  APPARATUS 

To  approximate  the  characteristics  of  an  infinite  conical  nozzle, 
a  nozzle  17.5  inches  in  length  with  an  8  inch  diameter  entrance  and  a 
1  inch  diameter  exit  was  designed  and  constructed.  It  subtended  a 
half-cone  angle  of  11.3  degrees  and  featured  an  entrance-exit  planar 
area  ratio  of  64.  This  ensured  that  stagnation  conditions  existed  at 
the  entrance  for  all  subsonic  flows  in  the  nozzle.  (See  Shaprio  [5], 
p.  614.) 

Figures  1  and  2  show  the  overall  arrangement  of  the  experimental 
apparatus  and  associated  electronic  equipment.     Figure  3  schematically 
depicts  the  test  section. 

Compressed  air  from  two  supply  tanks  passed  through  a  pressure 
regulating  valve  and  into  the  10  inch  diameter  stagnation  chamber 
through  a  baffling  arrangement.     The  stagnation  chamber  contained 
honeycomb  material  with  a  length/diameter  ratio  of  10  to  dampen  out 
any  disturbances  in  the  plenum  air.     A  fine  mesh  screen  was  also  installed 
to  hasten  turbulence  decay  prior  to  the  entrance  of  the  flow  into  the  nozzle. 
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Sinusoidal   pressure  disturbances  in  the  nozzle  exit  plane  were 
created  by  a  vibrating  piston  located  in  the  exhaust  chamber.     The 
piston  was  driven  by  a  MB  Electronics  Model   EA  1250  Vibramatic  Exciter 
powered  by  a  Model   2120  Amplifier  using  a  sinusoidal   signal   from  a 
Hewlett-Packard  Model   200CD  Wide  Range  Oscillator.     The  vibrator  was 
mounted  on  three  horizontal  tracks  so  that  the  zero  position  of  the 
piston  could  be  adjusted  for  optimum  operation. 

Dynamic  pressure  measurements  in  the  exit  plane  were  recorded  with 
a  Bruel  and  Kjaer  Model  4138  1/8  inch  Condenser  Microphone  with  Adaptor 
UA0036  and  Model   2618  1/4  inch  Preamplifier  connected  to  a  Model   2803 
Two-Channel   Power  Supply. 

Mass  flux  (pq)  measurements  were  performed  using  a  Thermo -Systems, 
Incorporated  Model   1050  Constant  Temperature  Anemometer  with  a  DISA     ■ 
Type  55F33  Hot  Wire  Probe.     The  hot  wire  anemometer  measured  both  steady 
state  (DC)  and  dynamic   (AC)  components  of  mass  flux.     The  anemometer  DC 
voltage  was  measured  on  a  Simpson  Model   2700  Digital   Voltmeter. 

The  anemometer  and  microphone  signals  were  passed  through  identical 
Krohn-Hite  Model   3750  Filters   (band-pass)  to  preserve  phase.     The  AC 
voltages  were  measured  by  two  Hewlett-Packard  Model   3400A  RMS  Voltmeters. 
The  filtered  dynamic  signals  were  displayed  on  a  Hewlett-Packard  Model 
132A  Dual   Beam  Oscilloscope  for  visual   verification  of  system  response. 

Phase  angles  between  pressure  and  mass-flux  fluctuations  were  meas- 
ured on  a  Dranetz  Model  305-PA-3002  Digital   Phase  Meter  using  the  output 
signals  from  the  oscilloscope  amplifiers.     Figure  4  shows  a  block  diagram 
of  the  electronic  measurement  equipment. 

Manometers  were  used  to  measure  stagnation  and  exit  static  pressures. 
Stagnation  temperature  was  measured  by  a  mercury-in-glass  thermometer. 
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B.  PROCEDURE 

Mass  flux  (pq)  calibration  of  the  hot  wire  anemometer  was  performed 
using  a  Thermo -Systems,  Incorporated  Model  1125  Calibrator.  The  method 

of  least  squares  was  used  to  fit  an  eighth  order  polynomial  in  the 

2     2 
quantity  (V  -  V  )  to  the  calibration  points.  (See  Computer  PROGRAM 

THREE.)  Microphone  calibration  was  performed  with  a  Bruel  and  Kjaer 

Piston  Phone:  Type  4220. 

Each  experimental  run  consisted  of  first  setting  a  steady  state 

flow  rate  by  maintaining  a  constant  pressure  in  the  stagnation  chamber. 

The  frequency  and  amplitude  of  the  pressure  disturbances  were  set  on 

the  oscillator.  Both  band-pass  filters  were  set  to  pass  only  a  narrow 

band  at  the  disturbance  frequency  with  gain  settings  of  20  db.  Data 

was  taken  over  the  range  of  frequencies  from  20  Hz  to  250  Hz. 

C.  DATA  REDUCTION 

Reduction  of  experimental  data  was  accomplished  on  a  digital  com- 
puter through  the  use  of  PROGRAM  FOUR. 

Steady  state  mass  flux  was  computed  by  using  the  anemometer  DC 
voltage  in  the  evaluation  of  the  hot  wire  polynomial.  The  magnitude  of 
mass  flux  perturbation  was  obtained  by  first  correcting  the  anemometer 
RMS  voltage  and  then  adding  it  to  the  DC  voltage  to  yield  a  peak  voltage. 
The  hot  wire  polynomial  was  then  evaluated  for  the  peak  voltage,  giving 
a  maximum  value  of  mass  flux.  The  steady  state  value  was  then  subtracted 
from  the  maximum  value  to  arrive  at  the  perturbation. 

Microphone  RMS  voltage  was  similarly  corrected  to  obtain  a  peak  value 
and  then  multiplied  by  the  microphone  sensitivity  to  yield  the  magnitude 
of  the  pressure  perturbation.  A  6  db.  narrow  band  insertion  loss  in  the 
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band-pass  filters  was  also  accounted  for  in  the  computations  of  dynamic 
quantities. 

Mach  number  was  computed  from  the  one  dimensional,  isentropic 
pressure  relationship.  The  velocity  q  used  in  the  computation  of 
Modified  Strouhal  number  S*  was  determined  by  dividing  the  mass  flux 
by  a  value  of  p  arrived  at  by  using  the  perfect  gas  law  and  the  isen- 
tropic pressure-density  relationship. 

V.  DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

Figures  5  through  14  show  a  comparison  of  the  experimentally  deter- 
mined amplitude  factors  and  phase  angles  with  values  predicted  by  the 
Modified  WKB  solution  for  Mach  numbers:  .084,  .117,  ,201,  .278,  and  .340, 
Both  amplitude  factor  and  phase  angle  are  plotted  against  modified 
Strouhal  number  S*.  The  smooth  curves  are  the  modified  WKB  solutions 
and  the  experimental  results  are  represented  by  crosses. 

For  low  frequency  sinusoidal  pressure  disturbances  with  correspond- 
ing small  values  of  S*,  the  phase  angle  is  -180  degrees  in  accordance 
with  the  modified  WKB  solution  which  asymptotically  approaches  that 
value.  An  increase  in  pressure  results  in  a  decrease  in  mass  flux  and 
physically  it  is  easy  to  visualize  an  effective  phase  angle  between  the 
rising  pressure  and  the  falling  mass  flux.  Mathematically,  this  effec- 
tive phase  angle  is  simply  180°  plus  the  value  of  the  actual  phase  angle. 
Thus,  for  low  frequencies  there  is  no  effective  lag  between  it  and  u. 

The  WKB  values  of  amplitude  factor  monotonically  increase  from  unity 
for  low  S*  to  the  asymptotic  limit  for  extremely  high  S*.  This  means 
that  the  magnitude  of  the  mass  flux  response  diminishes  as  the  frequency 
of  the  pressure  disturbances  increases. 
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Chiang  [2]  noted  that  the  quasi-steady  relationship  was  probably 
valid  for  low  values  of  S*  but  became  increasingly  inaccurate  for  larger 
values  of  S*. 

The  experimental  values  of  amplitude  factor  agree  quite  well  with 
those  predicted  for  all  Mach  numbers  tested.  The  experimental  values 
for  phase  angle  differed  from  those  predicted  but  exhibited  a  similar 
trend.  For  the  low  frequency  tested  (20  Hz)  all  experimental  phase 
angles  were  close  to  -180  degrees  while  the  WKB  values  were  slightly 
higher.  For  higher  frequencies,  the  experimental  phase  angles  became 
larger  than  the  WKB  predicted  values.  This  apparent  phase  shift  can 
be  attributed  to  effects  of  viscosity  in  the  flow.  Viscosity  serves 
to  dampen  the  flow  system  with  a  resultant  increase  in  effective  phase 
angle  over  the  undamped  system,  i.e.,  that  assumed  in  the  WKB  analysis. 

Physical  limitations  on  the  test  apparatus  prevented  investigations 
in  the  regions  above  M  =  .34  and  disturbance  frequencies  of  250  Hz. 
Electronic  noise  and  vibrations  generated  by  the  air  compressor  prevented 
its  being  run  during  experiments.  Consequently,  the  air  supply  in  the 
tanks  could  not  sustain  runs  over  a  frequency  range  at  higher  Mach  num- 
bers. The  other  problem  encountered  was  that  of  turbulence.  Early  in- 
vestigations revealed  that  turbulence  in  the  nozzle  was  prevalent  at 
frequencies  below  200  Hz.  Inclusion  of  the  fine  mesh  screen  in  the 
stagnation  chamber  helped  greatly  to  dissipate  turbulence  or  at  least  to 
break  it  up  into  high  frequency  turbulence.  Power  density  spectrums  of 
the  hot  wire  anemometer  signals,  obtained  by  taking  the  Fourier  Transform 
of  the  autocorrelation  function,  indicated  the  presence  of  turbulence 
with  frequencies  above  250  Hz.  This  high  frequency  turbulence  created 
too  high  a  signal  to  noise  level  for  meaningful  data  to  be  obtained 
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above  250  Hz.  A  PAR  Model  IOTA  Correlation  Function  Computer  and  a 
Model  102  Fourier  Analyzer  were  used  to  determine  the  power  density 
spectrums. 

It  is  concluded  that  the  modified  WKB  solution,  with  its  ideal 
flow  assumptions,  does  offer  a  reasonable  approximation  of  the  ampli- 
tude factor  for  pulsating  subsonic  flow  of  air  in  a  conical  nozzle, 
as  long  as  the  flow  originates  at  near-stagnation  conditions.  Pre- 
diction of  phase  angle  is  conservative  and  could  be  used  as  a  lower 
bound. 

It  is  recommended  that  further  investigations  be  pursued  at  higher 
Mach  numbers.  This  would  require  both  a  quiet  compressed  air  supply 
and  a  higher  power  vibrational  exciter.  It  is  felt  that  the  present 
location  of  the  air  compressor  within  the  laboratory  is  unsatisfactory 
and  chat  the  compressor  should  be  isolated  from  the  laboratory  by 
either  moving  it  outside  of  the  building  or  by  some  other  suitable 
means. 

It  is  also  recommended  that  an  analysis  be  undertaken  which  would 
account  for  the  effects  of  viscosity  upon  the  dynamics  of  a  pulsating 
flow  field  in  a  convergent  nozzle. 
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Figure  1.     TEST  APPARATUS 
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Figure  2.     ELECTRONIC  EQUIPMENT 
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FIGURE  5.   COMPARISON  OF  EXPERIMENTAL  AMPLITUDE  FACTOR  WITH 
•  MODIFIED  WKB  SOLUTION,  M  =  .084 
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FIGURE  6.  COMPARISON  OF  EXPERIMENTAL  PHASE  ANGLE  WITH 
MODIFIED  WKB  SOLUTION,  M  =  .084 
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FIGURE  7.   COMPARISON  OF  EXPERIMENTAL  AMPLITUDE  FACTOR 
WITH  MODIFIED  WKB  SOLUTION,  M  =  .117 
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FIGURE  8.   COMPARISON  OF  EXPERIMENTAL  PHASE  ANGLE 
WITH  MODIFIED  WKB  SOLUTION,  M  =  .117 
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FIGURE  9.     COMPARISON  OF  EXPERIMENTAL  AMPLITUDE  FACTOR 
WITH  MODIFIED  WKB  SOLUTION,  M  =   .201 
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FIGURE  11.   COMPARISON  OF  EXPERIMENTAL  AMPLITUDE  FACTOR 
WITH  MODIFIED  WKB  SOLUTION,  M  =  .278 
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FIGURE  12.   COMPARISON  OF  EXPERIMENTAL  PHASE  ANGLE 
WITH  MODIFIED  WKB  SOLUTION,  M  =  .278 
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